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I. INTRODUCTION 


The goal of the present research is to develop the capability to compute static 
aerodynamics on configurations such as guided and unguided projectiles for the Mach 
number range zero to three and angle of attack range zero to about twenty degrees. 
The Mach number and angle of attack range cover present and probable future 
design requirements for guntaunched weapons. 


Included in the present report is a detailed description of the computer 
program, along with several example cases and a FORTRAN listing of the program. 
For the derivation and discussion of the various theoretical methods used in the 
development of the prediction program, the reader is referred to Part I of this 
report (Reference |). However, for information purposes, the methods used to 
compute the particular force or moment component in the given Mach number 
region are listed in Figures | and 2. Figure ! gives the methods for the body alone 
and Figure 2 those for the tail (or canard) alone along with the interference effects. 
Most of the methods listed are standard in the literature with the exception of the 
empirical methods and the combined Newtonian Perturbation theories. Detailed 
discussion of this new theoretical method for calculating body wave drag can be 
found in Reference 2 and for wing wave drag in Reference 1. 
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il. PROGRAM DESCRIPTION 
A. Configuration Geometry 


The program is designed for four possible configurations: (1) wing alone, 
(2) body alone, (3) wing-body, and (4) canard-body-tail. Note that in present 
terminology. wing is interchangable with either a canard or tail. There are several 
different geometries which the wing or body may have as discussed below. 


1. Body 


The body may have a pointed or blunted nose. Blunt noses may have 
spherical caps or they may be truncated as shown in Figure 3. The program 
automatically accounts for this, however, if the correct inputs are given as decribed 
in the input section. In addition to being pointed or blunt, the nose may have up 
to two different ogive segments present. For example, on spin-stabilized projectiles 
there is norma!!: one ogive on the fuze and a different ogive between the fuze and 
shoulder. If the aerodynamics are desired in transonic flow, there is a minimum 
allowable nose length of 1.5 calibers due to the tcvle look-up procedure used there. 


The total body alone may end with the nose or it may continue with 
an afterbody. If an afterbody 1s present it is assumed to be cylindrical. Again due 
to the empirical estimation of aerodynamics in transonic flow, the afterbody must 
be less than ten calibers because this is the upper limit of the tables. Following the 
afterbody, a conical or ogival boattail may or may not be present. Instead of a 
boattail, a flare may be considered but the base drag must be disregarded because it 
is dernved for a bhoattail angle. Finally, the body alone may or may not have a 
rotating band present. 


2. Wing 


The wing is assumed to have one of two airfoil sections: a biconvex 
or modified double wedge. Both airfoil sections may have sharp or blunt leading 
and trailing edges. Also, the wing thickness to chord ratio and the slope of the 
airfoil section may vary all along the span. 


It may appear at first sight that assuming the airfoil section to be 
one of the two shapes above severly limits the program. This is not the case for 
projectiles and missiles. however, since the fin geometry normally is of a simple 
planform with no camber. Also, referring to Figure 4, the modified double wedge 
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FIGURE 4. TAIL (OR CANARD) INPUT GEOMETRY 
NOMENCLATURE FOR A MODIFIED 
BOUBLE WEBGE AMFOIL 


can cover all wedge airfoil shapes by adjusting the parameters CRIW and CR2W 
along with the angles GAW (1) through GAW (4). For example, if a double wedge 
airfoil is desired then it is sufficient to make the quantity [CRW - CRIW - CR2W) 
equal to zero and to make GAW (2) = GAW (3). 


B. Aerodynamics 


As mentioned previously, the various theoretical and empirical methods 
used to compute the static aerodynamics are listed in Figures ' and 2 and will not 
be discussed in this report as they are discussed fully in References | and 2. It is 
worthy of note that the overall guiding principle in the choice of methods listed in 
the above figures was to use techniques which could yield accuracies generally in 
the range +10% and which were inexpensive on the computer. It is believed that 
this goal was accomplished in that drag and normal force for most configurations 
can be obtained within the above accuracies and for a cost not exceeding $75.00 
for ten Mach numbers or angles of attack. The center of pressure is within a half 
caliber of experimental data for most configurations. 


C. Subroutines 


A brief description of each subroutine will aid the user in understanding 
the methodology and logic of the program. These descriptions follow. 


1. MAIN 


. This subroutine acts as a control for the entire program. It handles 
the calling of the various sutroutines that calculate the different lift and drag 
components. The calling of the various subroutines depend on the different options 
present in the program. These options are given in the input section of this report. 
The MAIN also handles the summation of the different components to obtain the 
total lift, drag, and pitching moment coefficients. 


2. AINTER 


A double interpolation routine used in the transonic Mach number 
range. 


3-8. ARCOSH, ARCOS, ARSECH, ARSIN, ARSINH, ARTANH 


These subroutines caiculate the hyperbolic arccosine, arccosine, 
hyperbolic arcsecant, arcsine, hyperbolic arcsine, and the hyperbolic arctangent, 


respectively. They are included because not all computers contain these particular 
functions in their library. 


9. BASEP 
Calculates the base drag for a body of revolution with or without a 
boattail throughout the Mach number range. It also includes an empirical estimate 
for the increase in base pressure drag due to the presence of fins. 


10. BASEPW 


Computes the drag due to trailing edge separation of a blunt trailing 


edge fin. 
11. BLUNT 
Derives the coordinates of a blunt nose tip on a body of revolution. 
12. CP3DW 


Determines the perturbation velocity at each point on a fin. These 
perturbation velocities are used by subroutine WING to compute pressures and forces 
en a fin. 


13-16. DISC1, DISC2, DISC3, DISC4 


These subroutines put in appropriate perturbation solutions to simulate 
discontinuities in body shape and curvature. 


17. DIST 


Computes the spanwise distribution of lifting pressure in subsonic 
flow. 


18-19. ELIPT1, ELIPT2 


These subroutines are used to evaluate the complete elliptic integrals 
of the first and second kind. 


20. FBINT 


Calculates the fin body and body fin interference. The method used is 
that of Nielsen and Kaattari.?) 


21. FOINT 


This is a simple linear interpolation routine. It is used in subroutines 
TRNCNA and SUBTRN. 


22-23. FDS, FOPS 


Theses subroutines are used to find the derivatives of a function at a 
given point. The five point Lagrange method is used. 


24. GCALC 


Computes the spanwise interpolation function for use in subsonic 
lifting surface theory. 


25. GEOM 


This subroutine reads input body coordinates and then computes the 
coordinates where the flow-field properties will be calculated. 


26. GEOM1 
_ Calculates the geometric properties of wings in subsonic flow. 


27. GUIDED 


This subroutine is the controlling subroutine for subsonic wing lift. It 
establishes boundary conditions and applies the Prandt] Glauert transformation. 


28. HCALC 


Computes the chordwise interpolation function used in subsonic lifting 
surface theory. 


29. HINT 


Computes the chordwise integrals to be used in the chordwise 
interpolation functions in subroutine HCALC. 


30. HYBRID 


Determines the body pressure using the hybrid perturbation theory of 
Van Dyke. For a discussion of this theory, see Reference (2). 


31-32. INTERP, INTERS 


These subroutines are used to interpolate for the value of a function 
at a given point. Five point Lagrange interpolation is used. 


33. LIFT 


Acts as an executive program responsible for the calculation of fin 
lift. It calls the interference subroutines and combines the lift of the isolated fins 
with the proper interference terms to obtain all fin lift components. 


34. MINVR 
Solves the matrix equation AX = B where A is a square coefficient 
matrix and B is a matrix of constant vectors. A! and |A/ are also available. 
Solution is by the Gauss-Jordan elimination method. 


35. NEWRAP 


Uses the Newton-Raphson method to solve for the mean skin-friction 
coefficient for a given Reynolds number and Mach number. 


36. NEWT 
Computes the pressure and static aerodynamics on the blunt portion 


of the nose using modified Newtonian theory. It also calculates the match point to 
combine Newtonian theory with perturbation theory. 
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37. NORMFO 


Solves for the normal force coefficients on the various components of 
the body in transonic flow using mostly empirical methods. 


38. PQORINT 
Integrates the upwash effect of one wing panel on another. 
39. RBAND 
Estimates the increase in drag due to the presence of a rotating band. 
40-41. REGONE, REGTWO 
These subroutines calculate the supersonic lift and center of pressure 
on a fin with a subsonic leading edge and a supersonic trailing edge by linear 
theory. 


42-46. REG1, REG2, REG3, REG4, REGS 


These subroutines calculate the supersonic lift and center of pressure 
on a fin with supersonic leading and trailing edges by linear theory. 


47. RK 


Finds the solution of differential equations using the tourth-order 
Runge-Kutta technique. 


48. SIMP 


Simpon’s rule is used to integrate surface pressures to find lorees and 
moments on a body. 


49. SIMPW & 


Intergrates wing pressure due to thickness to find wi. drag on a 
wing by using Simpson’s rule. 


Applies Mangler’s principle Value technique to obtain the solution to 
an improper intergral. 


51. SKINF 


Calculates the axial force coefficient due to skin friction on the ‘tody. 
62. SKINFW 


This subroutine computes the skin friction drag of a wing. The 
Reynold’s number is based on the mean geometric chord. 


53. SUBCNA 


Calculates the subsonic normal force coefficient and center of pressure 
for an isolated fin. It acts as a calling program for the subroutine that actually does 
the calculations. 


54. SUBTRN 


This subroutine calculates the transonic normal force coefficient and 
center of pressure for an isolated fin. The method usec can be found in the USAF 
Stability and Control DATCOM."4) The method used in the program is a sliht 
modification of that found in the DATCOM, but should be more accurate because 
the highest subsonic value of normal force derivative and lowest supersonic value are 
calculated by using lifting surface theory and linear theory, respectively. 


65. SUBXCP 


Calculates the subsonic center of pressure for an isolated fin. The 
subroutine uses values of sectional center of pressure as calculated in subroutine 
DIST. 


56. SUPCNA 


Computes the supersonic normal force coefficient for an isolated fir 
It acts as an executive program and it sets up fin geometry, determines what region 
a given point on the fin is in, and numerically adds all the lift increments in order 
to determine the normal force. The supersonic center of pressure for an isolated fin 
is also calculated in this subroutine. 
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57. TRANS 


Determines the wave drag of a boattail in transonic flow. It is also 
used to calculate the nose wave drag of tangent ogives in transonic flow. 


58. TRAPE 


Trapezoidal rule of integration used to determine the surface area and 
volume of the body alone. 


59. TRNCNA 


This subroutine is used, along with subroutine SUBTRN, to calculate 
the transonic normal force coefficient. It acts as an executive program as it calls the 
various subprograms necessary to determine a table of transonic normal forces for an 
isolated fin as a function of Mach number. 


60. WAVE 


Integrates the body alone pressures in order to compute the static 
aerodynamics. 


61. WING 
Calculates the pressures and forces on a fin due to thickness at 
supersonic .speeds. The airfoil thickness is assumed to be symmetrical about the 
x-axis. The method used is conical flow theory as modified in Reference 1. 
62. WTINT 
Finds the decrement in normal force derivative of the tail of a 


configuration due to downwash from the canards. The method used is that of 
Nielsen and Kaattari'?) mentioned earlier. 
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HI. INPUT 


The following is a list of the required inputs to the computer program 
described in this report. 


CARD TYPE | FORMAT (13) 
Variable Name Column Variable Description 
M (1-3) Number of cases to be run 


CARD TYPE !1 FORMAT (4F10.4,2F 15.12,215) 


Variable Name Column Variable Description 

AL (1-10) Angle of attack (Degrees) 

DIA (11-20) Reference diameter of body (Ft) 

HB (21-30) Mean height of the rotating band above the body 
surface (Calibers) 

AINF (31-40) Speed of sound (Ft/Sec) 

RHOINF (41-55) Density (Slugs/Ft?) 

AMUINF (56-70) Absolute viscosity (Ib-Sec/Ft?) 

IPRINT (71-75) Equal | if pressure coefficients are to be printed 
Equal 2 no pressure coefficients printed 

NTYPE (76-80) Equal | body alone aerodynamics calculated 


Equal 2 body-wing aerodynamics calculated 
Equal 3 body-wing-canard aerodynamics calculated 
Equal 4 wing or canard alone aerodynamics calculated 


CARD TYPE 11t FORMAT (6F5.3,15) 


Variable Name Column Variable Description 

XW (1-5) Distance of wing leading edge from nose tip 
(Calibers) see Figure 4 

DELTAW (6-10) Wing deflection angle (Degrees) 

DW (11-15) Diameter of body at wing root chord. If the 


diameter varies, an average of the body diameters 
at the leading and trailing edge should be used 
(Ft). 

XC (16-20) Distance of canard leading edge from nose tip 
(Calibers) 
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Variable Name 


DELTAC 
DC 


XCG 


MN 


Column 


(21-25) 
(26-30) 


(31-35) 


(36-40) 


Variable Description 


Canard deflection angle (Degrees) 

Diameter of body at canard root chord. If the 
diameter varies, an average of the body diameters 
at the leading and trailing edge should be used 
(Ft). 

Reference point for moments and center of 
pressure (measured in calibers from most forward 
point of nose). 

Number of Mach numbers to be computed 


CARD TYPE iV FORMAT (16F5.3) 


Variable Name 


AM(ARRAY) 


Column 


(1-80) 


Variable Description 


Mach numbers (Limited to 16) 


CARD TYPE V FORMAT (15F5.3,15) 


Variable Name 


GAW(1) 
GAW(2) 


GAW(3) 


GAW(4) 
CRW 
CTW 
BW 
CRIW 


CR2W 


Column 


(1-5) 
(6-10) 


(11-15) 


(16-20) 
(21-25) 
(26-30) 
(31-35) 
(36-40) 


(41-45) 


(46-50) 
(51-55) 
(56-60) 
(61-65) 


Variable Description 


Tail leading edge sweep angle (Degrees) 

Angle at which first line of sinks is swept back 
from Y-axis of tail (Degrees), see Figure 4. 

Angle at which second line of sinks is swept 
back from Y-axis of tail (Degrees) 

Tail trailing edge sweep angle (Degrees) 

Tail root chord (Ft) 

Tail tip chord (Ft) 

Span of isolated tail panels (Ft) 

Distance from _ tail leading edge to first 
discontinuity measured from root of tail parallel 
to freestream (Ft) 

Distance from tail trailing edge to first 
discontinuity ur:'ream from root chord parallel 
to freestream (Ft) 

Leading edge radius of tail at root chord (Ft) 
Leading edge radius of tail at tip chord (Ft) 

Tail thickness at root (Ft) 

Tail thickness at tip (Ft) 
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Variable Name 


XOC 


RTEW 


Column 


(66-70) 


(71-75) 
(76-80) 


Variable Description 


Distance of wing trailing edge from base (positive 
upstream of base and measured in root chord 
lengths) 

Trailing edge radius of tai! at root chord (Ft) 
Equal | double wedge or modified double wedge 
airfoil; Equal 2 biconvex airfoil 


CARD TYPE VI FORMAT (15F5.3,15) 


Variable Name 


GAC(1) 
GAC(2) 


GAC(3) 


GAC(4) 
CRC 
CTC 
BC 
CRIC 


CR2C 


RRC 
RTC 
TRC 
TTC 
XOCI 


RTEC 
IC 


Column 


(1-5) 
(6-10) 


(11-15) 


(16-20) 
(21-25) 
(26-30) 
(31-35) 
(36-40) 


(41-45) 


(46-50) 
(51-55) 
(56-60) 
(61-65) 
(66-70) 


(71-75) 
(76-80) 


Variable Description 


Canard Leading Edge Sweep Angle (Degrees) 
Angle at which first line of sinks is swept back 
from Y-axis of canard (Degrees), see Figure |. 
Angle at which second line of sinks is swept 
back from y-axis of canard (Degrees) 

Canard trailing edge sweep angle (Degrees) 

Canard root chord (Ft) 

Canard tip chord (Ft) 

Span of isolated canard panels (Ft) 

Distance from canard leading edge to first 
discontinuity measured from root of canard 
parallel to freestream (Ft) 

Distance from canard trailing edge to first 
discontinuity upstream from root chord parallel 
to freestream (Ft) 

Leading edge radius of canard at root chord (Ft) 
Leading edge radius of canard at tip chord (Ft) 
Canard thickness at root (Ft) 

Canard thickness at tip (Ft) 

Distance of canard trailing edge from base of 
projectile (positive upstream of base and 
measured in root chord lengths) 

Trailing edge radius of canard at root chi rd (Ft) 
“qual 1 double wedge or modified 

double wedge airfoil 

Equal 2 biconvex airfoil 


CARD TYPE Vil FORMAT (815,4F 10.5) 


Variable Name 


N 


NSHAPE 


Ni 


N2 


Column Variable Description 
(1-5) Total number of points 
to be read in to describe 
the body alone geometry (limit of 30) 
(6-10) Parameter used to describe 
the body geometry as defined below 
Pointed bodies 
NSHAPE=1 Nose only 
NSHAPE=2 Nose plus afterbody 
NSHAPE=3 Nose with a discontinuity (there may or may not 


be an afterbody present) 
NSHAPE=4 Nose plus afterbody plus boattail 
NSHAPE=5 Nose with discontinuity plus afterbody plus boattail 


If NSHAPF=3 or 5 at least five points must 
be read in along each of the nose sections, even if the nose section 
is a straight line. 


Blunted Bodies 


NSHAPE must be 3 or 5 

NSHAPE=3 NNIA=2 Blunted nose with a discontinuity so there are 
two nose sections present (no boattai!l present) 

NSHAPE=3 NNIA=!_ Blunted nose with no discontinuity (no 
boattail present) 

NSHAPE=5 NN!A=2 Blunted nose with a discontinuity so there are 
two nose sections present (boattail present) 

NSHAPE=5 NNIA=I Blunted nose with no discontinuity (boattail 
present) 

If NNIA=1, then NI=1) and N225 

If NNIA=2, then NI2S and N2>9 


(11-15) Number of points used to describe the first nose 
section 

(16-20) Number of points used to describe the second 
nose section plus the number of points used to 
describe the first nose section. 
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Variable Name 


N3 


RR 


Column 


(21-25) 


(26-30) 
(31-35) 


(36-40) 


(41-50) 


(51-60) 


(61-70) 


(71-80) 


Variable Description 


Equal | conical boattail 

Equal 2 ogival boattail (if ogival boattail is 
present at least five points must be used to 
describe the boattail section) 

Equal | pointed body 

Equal 2 blunted body 

Equal |, spherical cap on nose 

equal 2, truncated nose 

Equal |, no discontinuities present in nose 

Equal 2, discontinuity present in nose so nose 
appears to be made of two distinct sections 
Parameter used to describe mesh spacing. For 
blunted or spherically capped nose, C2=.05 and 
for pointed nose C2=.9 are nominal values 
Another parameter used to describe mesh spacing. 
For blunted or spherically capped nose, C4=1; 
for pointed nose, C4=20 are nominal values 
Constant which determines limiting body slope 
for a given mach number (.95 recommended) 
Radius of spherical cap or truncated nose 
(Calibers) 


CARD TYPE VII! FORMAT (2F 15.10) 


Variable Name 
x() - 
R(I) 


Column 


(1-15) 


(16-30) 


Variable Description 


Longitudinz] body coordinate measured from nose 
(calibers). If nose is blunt, X(1)=0 is at the end 
of the sp‘erical cap or at the truncated position. 
Body radius at given longitudinal station 
(calibers) 


There are as many Card Type VIII as arc needed to describe the body up to 


30 points. 


It should also be poioted out that if: 
NTYPE=1! (Body alone) Card Types V and VI are oniitted, 
NTYPE=2 (Wing-body) Card Type VI is left blank or contains all zeros, 
NTYPE=3 (Wing-body-canard) all Card Types contain data, 


NTYPE#4 (Wing only) Type VI is left blank and Card Types VII and VIII 
are omitted. 


If several cases are to be computed, Card Type I is input one time only, but 
Card Types II - VIII are given for each case as needed. 


Figure 3 shows a few typical examples of body geometry and associated 
parameters. 


At present, the maximum Mach number that can be input is three. This is due 
to limitations on tables present in the program. The upper limit on Mach number 
may be less than this if the local body slope becomes equal to 0.95 of the Mach 
angle based on the freestream Mach number. Mach numbers should be read in from 
highest to lowest. If a Mach number is desired between 1.05 and .°S, then 
Mach 1.05 must be included. This is because values of wave drag %n boattzils is 
assumed to vary linearly from 0 at M = .95 to an analytically calculated value at 
M = 1.05. Wave drag on fins is also assumed to vary linearly from 0 at M = 0.9 
to an analytically calculated value at M = 1.05. 


The shortest nose length that can be considered in the transonic flow regime is 
1.5 calibers. The longest afterbody length that can be considered in the transonic 
flow regime is 10 calibers. These limitations are due to limits in the tables that are 
internal to the program. 


Several sample input data sheets are shown in Figures 5 thru 8. An 
accompanying sketch with the geometric parameters describing each configuration is 
also. given. Figure 5 is a body alone, Figure 6 a wing-body, Figure 7 a 
canard-body-tail, and Figure 8 a wing alone. 
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FIGURE 5 — SAMPLE INPUT DATA FOR A TYPICAL BODY ALONE CONFIGURATION 
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Bata CAMO LATOUT 
fate Ot 
PO -OU - 10aRS 08 Don. 1-400 


SECTION B-B 


AR = 466,(t/c), =.1, (t/c), =.1,%*.5 


FIGURE S — SAMPLE NEPUT DATA FOR A TYPICAL WING ALONE CONFIGURATION 
( BICONVEX AIRFOM DESIGN) 


23 


IV. OUTPUT 


Before the aerodynamic forces and moments are given, the computer printout 
will list the freestream conditions and the wing-body geometries that were input. 
The force and moment output is then given in component form. Referring to the 
example computer output in Figure 9 (output corresponds to geometry input of 
Figure 7), the first table is the body alone drag which is broken down into 
skin-friction, base, wave, and protrusions (rotating band). The second and third 
tables are the tail and canard drags (if the canard is present), also broken down 
into the same drag components as above. The fourth table lists the normal force 
contributions from all geometry components: body alone, wing alone, wing-body, 
body-wing, canard alone, canard-body, body-canard, and canard-tail. The last table 
lists the total static aerodynamics of the entire configuration. These include drag, 
lift, pitching moment, and center of pressure. Also included are the secant slopes of 
normal force and pitching moment which for small angles of attack are the normal 
force and pitching moment coefficient derivatives. 
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FIGURE 9 — SAMPLE CUTPUT FOR CANARD-BODY-TAM CONFIGURATION OF FIGURE 7 


CASE NA, 1 


ANGLF OF ATTACK = el 0DEGS PFEFRENCE MTAMFTEP = .J250FT 


OFFFOFNCE CONODITIANS 


SPEED OF SNUND = 1116.°90 FY/SEC 
OFMITY e .0021769 SLUGS/FT3 
AMSOLUTE VISCOSITY © .000000376528 LA-SEC/FT2 


WING GEQMETRY(ONUALE WEOGE 02 MODIFIED DOUALF WEOGF AIPFOIL OFSIGN) 


SOAN= .ISOFT. 

ONOT CHNeDs .O96FT, 

TIP CH(PN= .O9KFT, 

LEADING ENGE SWFFPS20.000F%. 

FIPST LINE OF SINKS 215.000FGS. 
SFCONOD LINE OF SINKS=26.000FGS. 
TRAILING FOGE SWEEP 220.000FGS. 
FIPST CHOPD SEGMENT= .O33FT. 

PLEAD CHNOD SEGMFNT® JO3SFT. 

RONT THICKNESS =.0125FT, 

TIP THICKNESS? .O175FT. 

LFADING FOGE PANTUS AT RONT= ,QO2LFT. 
LEADING FOGE PANTUS AT TIO= .O021FY 
TRAILING FIGE PLUNTNESS= .0062FT 
OEFLECTION ANGLE 0.00DEGS. 


CANARD GEOMETRY(DOURLE WEDGE 0° MODIFIED DOURLE WENGE ATRFOTL DESTGN) 


SPAN= .LATFT. 

ROOT CHNRD= .125FT, 

TI? CHOWN: .O79FT. 

LEADING EDGE SWEFPx 30. 000FS. 

FIPST LINE OF SINKS 216.00NEGS. 
SFCOND LINE OF SINKS=16.000FS. 
TRAILING FOGE SWEFP=s 0.000F6S. 
FIPST CHNPN SEGMENTS .062FT, 

PEAR CHIPD SEGMENT® .O42FT. 

PNAT THICKNESS: .0125FT, 

TIP THICKNESS= .OL25FT. 

LEADING EDGE PADTIUS AT RONTe ,O021F 7. 
LEADING EDGE PADIUS AT TIP= ,O021FT 
TRATLIMG EDGE ALUNTNESS= .0062FT 
OFFLECTION ANGLE 0.000EGS. 


BODY CONPOTNATES 


x Rr 

0.6000 20936 

0 2 UG 01718 

24987 02616 

o THA 3028 

©9973 3556 
12467 3999 
1.3317 06138 
1. & 167 46266 
1.25016 ob367 
1.5867 74498 
1.6716 «46603 
1.767 +4698 
te AbLt6 04786 
1.9267 ob AGS 
2.0116 04937 
2.0%7 5000 
& 6700 «5000 
5.3400 23949 
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wach NO. 


2.060 
2.400 
2.06¢ 
1.600 
1.200 
1.05¢ 
1.000 
950 
4900 
«700 
«500 


watw NO. 


>.a0l 
2.600 
2.000 
1.40¢ 
1.200 
1,050 
1.000 
0980 
2900 
2700 
5tC 


macH NO, 


7.50C 
2.600 
?.occ 
1.6090 
1.200 
1.00 
1.900 
3 9F0 
©900 
700 
500 


macy NO, 


2.0000 
2.6009 
2.0009 
1.6000 
1.2009 
1.0560 
1.0000 
©9500 
9000 
7000 
«5080 


SMIN 


SKIN 


ANDY SLONE 


wace Wn, 


20 600 
20600 
2.600 
1.660 
12200 
1.050 
1-600 
° 960 
2900 
0700 
$00 


0087 
10056 
004A 
20040 
.0032 
-003t 
20031 
2003! 
.0otr 
0027 
0026 


paic 


0031 
2036 
037 
«O41 
2 Obs 
2066 
2046 
20a7 
0067 
0048 


TION 


8 
6 
9 
3 
a 
2 
6 
1 
5 
a 


20095 


Fore 


2006 
2006 
«007 
2007 
2006 
006 
20046 
2096 
2005 
2008 
20046 


Fetc 


003 
+003 
+003 
003 
+003 


TION 


5 
a 
1 
iC, 
q 
6 
LJ 
3 
1 
q 
3 


TION 


3 
5 
6 
? 
6 


20035 


2008 
008 
0003 
2002 
2003 


WING ALONE 


en 


04119 
ehhO? 
5016 
587° 
7571 
24560 
FOAL 
5050 
« 3836 
01935 
01067 


i) 
3 
2 
? 
1 


09017 
20021 
+0026 
20035 
20056 
2 106s 
0039 
20027 
«0021 
200466 
20042 


BNOY AYTAL FORCE CON 'RIAUTIONS 


PASE SRESSIIOF 


O71 
-0%2 
1193 
08555 
0121 
01013 
01736 
01565 
21322 
01102 
01019 


WING ANTAL FORCE CONT@TAUTT ONS 


PASF SPF SSUPF 


20041 
00105 
0137 
o01A8 
28266 
20279 
20276 
rr 
00207 
20379 
20175 


PRES SUPE 


of ele 
016493 


ORE SSUCE 


09% 
21055 
1160 
138 
22368 
22977 
22233 
01689 
oO Tey 
0.0000 
6.0009 


CANE@O ANTAL FORCE CONTOTAUTT ONG 


PASE OPE SSPE 


20039 
0050 
00165 
-0c84 
2126 
oO 1T8 
oCitt 
00118 
0099 
2004S 
of C86 


ORF SSUCF 


oO SAR 
20411 
20456 
00537 
00722 
01009 
0780 
0509 
0250 
9.0000 
0.0090 


NOPPAL FORCE CONTOTRUTIINS 


Canaen alone WING-RONY 


200cer 
«0020 
00012 
0015 
20020 
0013 
De 
00013 
001% 
20026 
0013 


0022 
20025 
0031 
00962 
20069 
20058 
20068 
0032 
20926 
«0056 
0051 


ANN WING cana@n-nony 


©6000 
90001 
0001 
20002 
0005 
20015 
0015 
20010 
2000" 
00017 
00016 


00013 
o OO1F 
2001" 
2802" 
oO 
21990 
2097? 
20071 
00971 
2007? 
00971 


TOVAL STATIC AE@PONYNAMICSIFOOCE sal Onn) 


cL ce 


CwaL 


$.233 
$.368 
5.567 
5.937 
7.298 
7.037 
6.932 
$.00¢ 
5.010 
7.253 
6.623 


cHaL 


“13.665 
“16.023 
716. 6Ru 
~1€,017 
22.936 
-216506 
20.231 
-16.710 
17,676 
-23.379 
“21.°53 


eenteu 


0.0000 
0.0¢000 
6.ceoe 
¢.caee 
0.0906 
6.0000 
0.0000 
o.coce 
6.0000 
¢.0000 
o.ceee 


TATAL 


e387 
212728 
eo 1 THA 
2 1638 
0 270t 
03503 
02573 
01799 
©1012 
20232 
oC23e 


TAT aL 


oChtO 
oCQ9R 
ocses 
20662 
oOAPd 
0116" 
0918 
oCKEY 
o0SA1 
C112 
0126 


5007 
-toce 
cores 
- 00f% 
008 
6023 
2 0016 
0023 
00013 
2 0016 
0013 


SIONS 


PADY-CANBOD = CANAPD-WING 


- 20002 
2003 
~ 20005 
-.ccor 
-.0017 
20009 
-.0c0r 
~ 0008 
~-CCO% 
~ooeie 
-.0t09 


¥C#70 


2.607% 
206126 
2.6671 
207822 
3.1327 
3.0569 
2.9285 
20*337 
2.3337 
3.2231 
3.2028 


TOTAL 


27522 
o27et 
«3093 
1876 
3987 
4070 
3397 
«7599 
020638 
21590 
otSi6 


YOvaAL 


20091 
= 609% 
00097 
0106 
00127 
20123 
0121 
00102 
20098 
0127 
20119 


V. COMPARISON WITH EXPERIMENT 


Three cases are considered to show, first of all, the general accuracy of the 
method when compared with experiment and, second, how the program can be used 
to obtain engineering estimates of aerodynamics for configurations which do not 
exactly fit into one of the four categories listed previously. The three configurations 
are a body alone, a wing-body in which the wing is mounted on a strake, and a 
canard-body-tail in which the tail does not have streamwise tips. 


The first of these configurations is the 5”/54 Rocket Assisted Projectile (body 
alone) for which the configuration geometry and input data are given in Figure 5 
and the aerodynamics in Figure 10. This particular spin stabilized projectile has a 
nose length of about 2.5 calibers and a_ boattail length of 0.5 caliber. The 
theoretical drag coefficient is in very good agreement with experiment throughout 
the Mach number range. Fair agreement is obtained for normal force coefficient 
derivative and center of pressure. The normal force coefficient derivative is generally 
low in the lower supersonic speed range and approaches the experimental data at 
moderate supersonic Mach numbers. 


A_ wing-body configuration is shown in Figure 6 and the _ corresponding 
aerodynamics in I1A and 11B. Note that in Figure 6. the tails are mounted on 
strakes which raises the question, “What does one use for the base diameter and 
wing planform?” The base diameter mainly determines the afterbody drag and the 
wing planform the wing lift. Since both of these quantities are a direct function of 
the base area and wing area respectively, it seems reasonable to compute these areas 
and then define the base diameter and wing planform from them. Thus the base 
area is computed inch:4'1g the strakes and then an equivalent base diameter defined. 
Next, the total wing «cco including the strake planform area, is computed and an 
equivalent wing obtained by adding this additional area to the chord and span. 
Although this does not change the configuration lift appreciably when based on 
wing area, it does change the lift considerably here because the wing lift is based 
on the body cross-sectional area. Using the above geometry modifications, the 
aerodynamics were computed and compared with experimental data in Figures 111A 
and 11B. The drag and center of pressure are shown in 1!A and the normal force 
coefficient denvative in 11B. Excellent agreement with experiment is obtained for 
normal force and center of pressure. The theoretical drag is about ten percent high 
at transonic Mach numbers but according to Reference 5, the blockage of the test 
model in the wind tunnel was too high. Normally when the wind tunnel model is 
too large (too much blockage), drag values measured in the transonic Mach range 
fall off and it appears that this may have happened in this case to account for 
some of the above discrepency. 
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FIGURE 10 


Comperison Theory and Test Data for 5°/54 Rap Projectile 
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FIGURE 11(A) 


Drag and Center of Pressure for a 
Typical Missile Configuration; AR = 4.5 
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FIGURE 11(B) 


Normal Force Coefficient Derivative for a Typical Missile Configuration; AR=4.5 
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The third geometry considered is a very general canard-body-tail shown 
schematically in Figure 12A along with the normal force and center of pressure. 
Note that the tail does not have streamwise tips so the total wing area is again 
computed and an equivalent span calculated based on the chord and wing area. This 
gives reasonable values for normal force and center of pressure (Figure 12A), 
although it appears the tail alone lift is about ten percent too high causing a 
rearward shift in the center of pressure. The drag of the body alone, along with the 
canard and tail components, is shown in Figure 12B. The body alone drag agrees 
well with experiment in subsonic and supersonic flow but is unacceptable in 
transonic flow. This is because of the sixty percent blunt nose which the empirica! 
transonic drag methodology does not account for. The wing drag is also high in 
transonic flow, but this is caused by the increase in body base pressure due to the 
presence of tail surfaces. This increase in drag is included in the curves at the 
bottom of Figure 12B. The empirical estimates of the body base pressure change 
due to fins is much higher for this case than the data suggest. However, the total 
configuration drag agrees with experiment within the accuracy bounds previously set 
forth, except in transonic flow. 
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FIGURE 12(A) 


Normal Force and Center of Pressure of 
A Missile Configuration; Rx, a4, R. =2, M_ =1.6 
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FIGURE 12(B) 


Drag of a Missile Configuration and Its Components 


33 


REFERENCES 


I. 


a 


Moore, F. G., Aerodynamics of Guided and Unguided Weapons: Part I — Theory 
and Application, NWL Technical Report TR-3018, December 1973. 


Moore, F. G., Body Alone Aerodynamics of Guided and Unguided Projectiles at 
Subsonic, Transonic, and Suppersonic Mach Numbers, NWL Technical Report 
TR-2796, November 1972. 


Nielsen, J. N. and W. C. Pitts, Wing-Body Interference at Supersonic Speeds 
With an Application to Combinations With Retangular Wings, NACA Technical 
Note TN-2677, 1952. 


Douglass Aircraft Co., Inc., USAF Stability and Control DATCOM, Revisions by 
Wright Patterson Air Force Base, July 1963, 2 Vols. 


Touch, L. M., Transonic Wall Interference Effects on Bodies of Revolution, 
AIAA Paper No. 72-1008. 


Glossary 


Aspect ratio of canard 

Aspect ratio of tail 

zero lift drag coefficient 

Lift coefficient 

Pitching moment coefficient 

Pitching moment coefficient derivative 
Normal force coefficient 

Normal force coefficient derivative 
Freestream Mach number 

Leading edge radius of tail or canard (ft) 
Trailing edge radius of wing (ft) 
Thickness to chord ratio of wing at root 
Thickness to chord ratio of wing at tip 
Coordinates with x along body and y out right wing 


Center of pressure measured in calibers from nose tip unless otherwise 
specified 


Angle of attack 
Canard deflection angle 


Wing taper ratio 
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